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charge and spin excitations; and c) a non- Fermi liquid state with incoherent charge 
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and asymptotically decoupled charge and spin excitations. The non-Fermi liquid 
phases occur for a wide range of parameters. 
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I. INTRODUCTION 



Various anomalies in the normal state of the high copper oxides have motivated 
the current theoretical efforts to search for models which exhibit metallic non-Fermi liquid 
states. The mechanism for the breakdown of Fermi liquid theory in one dimension is 
well understood. For weakly interacting one dimensional fermion systems, the renormaliza- 
tion group leads to the g-ology classification of spatially homogeneous metallic states. The 
possible states are Luttinger liquids or those with divergent CDW, SDW or superconducting 
correlation functions. Bethe Ansatz, bosonization and exact diagonalization methods 
have been used to show that this classification persists for strong interactions provided that 
phase separation does not intervene. |^ 

In higher than one dimensions, perturbative renormalization group analysis has shown 
that, Fermi liquid theory does describe weakly interacting fermion systems with a regular 
density of states. Therefore, the mechanism for the breakdown of Fermi liquid theory is 
necessarily non-perturbative in the interaction strength. Most of the previous approaches to 
this problem address under what circumstances the residual interactions between the quasi- 
particles can acquire singularities. Here we take an alternative approach. We analyze the 
competition between local correlation and itinerant effects without assuming quasiparticles 
to begin with. The effect of interactions is treated non-perturbatively. This is carried out 
in a controlled fashion, by studying this problem in the limit of infinite dimensionality. 

In the limit of infinite dimensions, all correlation functions of a lattice fermion model can 
be reconstructed from the knowledge of the local correlators, which can in turn be determined 
from those of an associated impurity model. P] The impurity model describes some local 
degrees of freedom interacting with a self-consistent electron bath. This limit corresponds 
to a dynamical mean field theory, within which non-trivial phases can be studied. The fact 
that only the low frequency behavior of the correlators of the bath governs the low frequency 
asymptotics of the correlators of the impurity model, suggests that a complete classification 
of the translationally invariant phases of fermions in the limit of large dimensions could 
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be within reach. The Hubbard model with a Lorentzian density of states maps onto an 
Anderson model with an infinite bandwidth, proving that this model in large dimensions is 
a Fermi liquid. |p The Hubbard model in large dimensions with a bounded density of states, 
at half filling and for large U, maps onto an Anderson impurity embedded in an insulating 
medium |TD|,|TT],0], demonstrating that at large U a frustrated Hubbard model is a Mott 



insulator. The characterization of the state that results from doping this Mott insulator is 
still an open problem. 

In this paper, we construct and study various metallic non-Fermi liquid states. These 
phases are realized in an extended (two band) Hubbard model which contains the most 
general local density- density and spin-spin interactions and a Lorentzian density of states. 
We stress however that our low energy analysis of the non- Fermi liquid phases presented here 
is more general than the particular extended Hubbard model that we studied. We envision 
that many realistic Hamiltonians will be described, at low energies, by the impurity actions 
that we study in this paper. The spin and charge couplings we discuss will be generated 
dynamically in the process of renormalization. 

The impurity model associated with the extended Hubbard model is a generalized An- 
derson model. It will be shown that the fact that the impurity model is associated with 
a lattice model forces the impurity to be near criticality. Qualitatively, this can be under- 
stood very simply in the context of a weakly doped large U Hubbard model. In this case 
the correlated site describes the local physics of the one band Hubbard model and the bath 
describe the effects of the electrons near that site. The self consistency condition restores 
translational invariance and forces the correlated level to be near the chemical potential. 
In a two band model describing a heavy (correlated) level and a light conduction band the 
heavy level participates in the charge transport only when it is situated near the chemical 
potential of the conduction band. The low energy description of the local physics of the one 
and the two band models is very similar. The role of the light conduction band in the two 
band model is played by the incoherent part of the one particle Green's function in the one 
band model. 



To study the low energy behavior of the correlators of the impurity model near the critical 
point, we set up a renormalization group analysis for this problem. Our strategy is to set 
up an expansion in terms of the hopping amplitudes between the local configurations of the 
impurity problem. This is an extension of the classic work of Haldane. [|13| The resulting 
expression of the partition function resembles that of a classical Coulomb gas. The hopping 
amplitudes correspond to fugacities of the Coulomb gas representation. The fugacities are 
proportional to the amplitude for making a transition between two different atomic states. 
When the amplitude for making transitions between states of different charge grows as we go 
towards longer time scales the system is a Fermi liquid. When this amplitude renormalizes 
to zero, quantum coherence is destroyed and Fermi liquid theory breaks down. We cast the 
breakdown of Fermi liquid theory in the framework of the macroscopic quantum tunneling 
(MQT) problem. [0] The transitions between Fermi liquid and non-Fermi liquid phases are 
extensions of the well known localization transitions in the MQT problem with one essential 
difference (responsible for a richer phase diagram) that we deal with a special three level 
system describing the local spin and charge degrees of freedom instead of the canonical two 
level system studied in the MQT literature. 

The possible phases are determined from the possible fixed points of the RG. We carry 
out the analysis of a more general SU(N) impurity model. The case N=2 corresponds 
to the three level system (empty site and SU(2) spin). We identify, besides the strong 
coupling phase in which all excitations are coherent and the weak coupling phase in which 
all excitations are incoherent, a new class of metallic state of matter in which spin excitations 
are coherent while charge excitations are incoherent. We call this an intermediate coupling 
phase. Both the intermediate coupling and weak coupling phases correspond to incoherent 
metallic states with vanishing quasiparticle residue, self-similar local correlation functions, 
and asymptotically decoupled charge and spin excitations. The low energy behaviors are 
not described by Fermi hquid theory. 

The N=l case reduces to the well studied resonant level model. flTRM It is well known 



that orthogonality effects are more pronounced when there are bound states which can 



violate the unitarity bounds of the couphng constants or in the presence of many channels 



of conduction electrons. This later point is emphasized by Ruckenstein and Varma in 
the search of a mechanism for Marginal Fermi liquid states in realistic models. Along this 
line, Giamarchi et al. have recently suggested that an impurity model containing a bath 
with several channels of screening electrons is closely related to the copper oxides. [0 More 
recently Perakis et al. suggested that the presence of several channels of screening electrons 
can give rise to non- Fermi liquid behavior similar to the one found in the multichannel Kondo 
problem. The addition of extra channels of conduction electrons ||20| in our formalism 



amounts to a change in the initial conditions of the renormalization group flow. 



The paper is organized as follows. In Section P A| , we introduce the extended Hubbard 
model and derive, in the limit of infinite dimensions, the associated impurity problem in a 
self-consistent medium. This impurity problem is a generalized asymmetric Anderson model. 
The qualitative physics associated with such an impurity problem is discussed in Section 
p.1 B| . The justification for the associated impurity problem being near criticality is given in 



Section [11 (J| . In Section [111 A| , we construct the atomic expansion for the partition function 



of the impurity problem. The resultant partition function is a summation over histories 
of the local degrees of freedom. It is shown in Section [111 B| that, this partition function 



corresponds to a particular realization of a class of one dimensional classical spin models 
with long range interactions considered by Cardy, extended to account for a symmetry 
breaking field. Using Cardy's procedure, we derive in Section |V| the renormalization group 
equations, and analyze the resulting universality classes and crossover behavior. In section 
we establish the nature of the Fermi liquid and the non- Fermi liquid phases of the extended 
Hubbard model. A summary of the phase diagram is given in Section |VC| . In Section |VD| , 
we discuss the relevance of our results to more realistic densities of states. In section [V^ 
we conclude with a discussion of the relevance of our results to realistic finite dimensional 
systems. The paper includes three appendices. Appendix A summarizes the bosonization 
procedure. Appendix B gives the details of the derivation of the renormalization group 
equations. Finally, in Appendix C, we derive the correlation functions in the weak coupling 



mixed valence regime. A brief account of the renormalization group flow and its implication 
for the extended Hubbard model in infinite dimensions has been given in Ref. p^ . 



II. THE EXTENDED HUBBARD MODEL IN INFINITE DIMENSIONS 

A. The Extended Hubbard Model in Infinite Dimensions 
and its Associated Impurity Model 

We study the following extended Hubbard model 

N ^ U 

ij (7=1 i (7=1 i a^cr' 

^ V V 

+ K4aCi(7 + h.C.) + ^Y.Y1 4(7CiAa'di<7' + ^Y.Y. 4aCia'dl,dia (1) 

i a=l i a, (7' i G^cr' 

It includes two species of electrons. The uncorrelated conduction c electrons have the hop- 
ping matrix tij or, equivalently, the dispersion e^. The localized d electrons have an energy 
level e^, and an on-site Coulomb interaction IJ which is taken to be infinity. These two 
species of electrons are coupled, at every site, through a hybridization t, a density density 
interaction and a spin exchange interaction For generality, we work with arbitrary 
spin degeneracy A^, which will also allow us to make contact with results from large 
expansion. Finally, the chemical potential is /x. 

In the spinless case (N=l) this model is reduced to a generalized Falicov-Kimball model 
which includes, besides the usual interaction Vi, a hybridization t term. In the limit of 
vanishing hybridization, the Falicov-Kimball model is exactly soluble in infinite dimensions. 
p3| , |2^ When the density of states is a Lorentzian p5|J^ all the correlation functions can 
be computed analytically and the origin of the breakdown of Fermi liquid theory in the 



model can be understood in simple physical terms. |^ There is a range of electron densities 
where the chemical potential coincides with the position of the renormalized d electron level. 
We referred to that situation as the pinning of the chemical potential. We notice however 
that the Falicov-Kimball model should be understood as the zero hybridization limit of the 



model (|T]). [^] In this paper we show that there are regions of parameter space in which 
the hybridization dynamically renormalizes to zero at low energies and establish that the 
Falicov-Kimball fixed point in the pinned region has a finite basin of attraction. 

The spin | case (N=2) is a version of the extended Hubbard model first studied in 
the context of the mixed valence problem. It has some similarities with the extended 
Hubbard model relevant to the high T^, systems proposed by Emery and Varma, Schmitt- 
Rink, and Abrahams. |^ 

The infinite dimension limit of a lattice fermion model is defined through an appropriate 
scaling of the hopping term, which in our case is that of the conduction electrons. It has 
been shown that, one can naturally associate with each lattice model in infinite dimensions 
an impurity model. This impurity model is obtained by integrating out all the degrees of 
freedom except for those living at the origin. 

Following this procedure we obtain the single site effective action of our extended Hub- 
bard model: 
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where tf)'^ = (c^,(i^). The Weiss effective field is a 2 x 2 matrix. It describes the effect 
of all the Feynman trajectories in which the electron at site zero leaves the origin, explores 
the lattice, and returns to the origin. Translational invariance demands that all the local 
correlation functions of the lattice model are the same as the correlation functions of the 
impurity model leading to the following self-consistency equations: 
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Here the non-interacting density of states, NQ{e) = J2k^{^ ~ ^fe) where the conduction 



electron dispersion, is the only place where the nature of the lattice enters. S*'"^ is a 
functional of Go defined by the self-energy matrix of the single-site impurity model defined 
in Eq. (^. After solving the system of equations (2) and (3) for Go we can evaluate the 
self-energy of the lattice model S = T,imp{Go)- 

For a Lorentzian density of states, Go can be solved explicitly 
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where F is the width of the Lorentzian. A more general density of states, however does 
not alter the qualitative results derived here, provided that the system is metallic and the 
hybridization t is smaller than the bandwidth of the conduction electrons. This will be 
discussed further in Section |V D[ 

Our single site action describes an impurity model : 
r0 
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Simp = drY^ 4<7(^r + /U - ek)cka + dl{dr + fi- ed)da + t{dlc„ + h.c.) 
+ ^ 4ca4'^<x' + ^ II c\c^>d\,d^ dld^d\,d^> 



(5) 



where an electron bath, with dispersion e^, has been introduced to generate the Weiss field 
G„ ^ given in Eq. (^. This is the generalized Anderson impurity model. 

B. Qualitative Physics of the Generalized Anderson Model 

The Hamiltonian of the generalized Anderson impurity problem corresponding to the 
action in Eq. @ is given by: 



^^imp = Y ^kcl^'^ka + E^dld^ + Y ^(4ca + h.c] 



fccr 



+^ cicadlda' + ^ II clc^'dld^ did^dldo 



(6) 



where the d electron level E'^ is given by 



E"d = ^d-^^ 



(7) 



and the bath electron dispersion Ek = — At- 

In this subsection, we give a quahtative discussion of the physics associated with this 
model. We separate the discussion for the spinless and spinful cases. 

In the spinless case (N=l) this model is reduced to the resonant level model. It was pre- 
viously established that, the resonant level model can be mapped onto an anisotropic 
spin I Kondo problem through the following: a) the hybridization t is identified with the 
transverse Kondo exchange b) the density- density interaction Vi, when written in a 
particle-hole symmetric form, is related to the longitudinal component of the Kondo ex- 
change J|| through identifying the corresponding phase shifts: (1 — 25/7r)^ — > 2(1 — 25||/7r)^, 
where 5 = tan~^ {it poVi/ 2) and 6\\ = tan~^{TxpoJ\\/^), with po the conduction electron den- 
sity of states at the Fermi level; and c) the d electron level E'^ is mapped onto the magnetic 
field {—H) applied on the local moment. Such a mapping allows a parallel discussion of the 
qualitative physics of the resonant level model and that of the Kondo problem. 

We recall that, at zero magnetic field, the Kondo problem has weak coupling fixed points 
when the exchange coupling is ferromagnetic: the transverse Kondo exchange is renormalized 
to zero. For an antiferromagnetic exchange coupling, it flows to a strong coupling fixed 
point. PO,PT| The transition can be viewed as an unbinding of defects {spin flips) in the 



path integral representation fo the partition function to be derived in section |T|, as in the 



Kosterlitz-Thouless transition of the 2D XY model. |^ In complete analogy, when the 
d level is located at the chemical potential, the resonant level model has a weak coupling 
fixed line for a range of attractive interactions given by 26 /n < —{V2 — 1). ||l5| , p!6| The 
hybridization becomes relevant for 26 /n > — (-\/2 — 1), at which a strong coupling phase 
emerges at low energies. The zero-temperature transition between these phases involves the 
unbinding of defects describing transitions between states with different charge {charge flips) 
in the path integral representation of the partition function to be described in section 

Applying a large magnetic field at the impurity in the Kondo problem polarizes the 
impurity spin, cuts off the infrared singularities, and therefore destroys the Kondo process. 
Similarly in the resonant level model, when the d level is far away from the chemical potential, 
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the d orbital is either fully occupied or empty, and the resonant level model becomes a 
potential scattering problem. Interesting physics occurs in the mixed valence regime, when 
the renormalized d-level, denoted as E^, is close to zero. 

The nature of these fixed points can also be discussed in parallel. In the ferromagnetic 
Kondo problem in zero magnetic field, the spectrum is composed of conduction electrons and 
a local spin mode, which are asymptotically decoupled at low energies. There is no energy 
scale in the problem; Correlation functions have long time algebraic behavior and do not 
have the Fermi liquid form. Similarly the spectrum of the resonant level model in the weak 
coupling regime, when the d-level is located at the Fermi energy of the conduction electrons, 
is composed of the local d-electron and the conduction electrons which are asymptotically 
decoupled. The d electron spectral function, as well as local correlation functions, exhibit 
power law behavior in frequency with interaction-dependent exponents, due to the finite 
value of the renormalized density-density interaction V*. The divergence occurs at zero 
frequency when — 0. In both cases, the low energy excitations are incoherent, and can 
not be analytically continued to those of non-interacting electron systems. 

In the antiferromagnetic Kondo problem, the local spin mode is quenched below the 
Kondo energy scale through the formation of quasiparticles, which correspond to co- 
herent spin excitations. In the resonant level model at the strong coupling regime, the 
renormalized hybridization t* between the local d and conduction electrons is finite. This 
renormalized hybridization sets the energy scale, below which the local d electron retains its 
free-electron character and combines with the conduction electrons to form quasiparticles. 
These quasiparticles describe coherent charge excitations. In both cases, the low energy 
behavior is described by a Fermi liquid. 

The above discussion indicates that, whether or not the hybridization is relevant deter- 
mines the nature of the low energy behavior. This point can be further seen by noticing 
that the wave function renormalization factor associated with the d electron is equivalent 
to the rcnormahzation of the hybridization: therefore, whether or not the the renormahzed 
hybridization is non-zero determines the existence or vanishing of the quasiparticle residue. 
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This motivates us to study the flow of the hybridization, along with other couphngs. 

For spinful case (A^ > 2) the generahzed Anderson model has low energy excitations 
in both the charge and spin channel. We will show, using the RG flow, that the interplay 
between charge and spin dynamics lead to three kinds of mixed valence states, characterized 
by different ways the hybridization and the exchange coupling are renormalized. 

Strong coupling mixed valence fixed point. Both the hybridization and exchange coupling 
are relevant. As in the antiferromagnetic Kondo problem and in the strong coupling regime 
of the resonant level model, the spectral functions at low energies are described in terms of 
a resonance around the chemical potential. The width of the resonance sets the scale for 
the low energy coherent quasiparticle excitations. This is in the same universality class as 
the usual strong coupling phase of the Anderson model. 

Weak coupling mixed valence fixed points. Both the hybridization and the exchange 
coupling are irrelevant. The system is self-similar; Correlation functions in the charge and 
spin channels have algebraic behavior similar to those of the ferromagnetic Kondo problem 
and of the weak coupling regime of the resonant level model respectively. Both the charge 
and spin excitations are incoherent. 

Intermediate coupling mixed valence fixed points. Here, the hybridization is irrelevant, 
while the exchange coupling is relevant. The system is self-similar in the charge channel, and 
exhibits a flnite energy scale in the spin channel. The correlation functions in the charge 
channel have an algebraic behavior similar to those of the weak coupling regime of the 
resonant level model, while the correlation functions in the spin channel are characterized 
by a resonance in analogy to those of the antiferromagnetic Kondo problem. Therefore, 
the low energy charge excitations are incoherent while the low energy spin excitations are 
coherent. 

The intermediate mixed valence phase occurs because the two kinds of defects (spin flips 
and charge flips) can unbind at different stages. |^ We will establish that, due to the special 
form of the coupling between the spin and charge dynamics, a phase with coherent charge 
excitations and incoherent spin excitations can not exist. 
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C. Pinning of the Local States at Chemical Potential 

In this subsection we show that, the fact that the impurity model is the result of a 
mapping from a lattice model forces the impurity to be near criticality. 

The local Green's functions of a lattice fermion model in infinite dimensions are given by 
those of the impurity problem. Therefore, the density of d electrons plus the local density 
of c electrons in the ground state of the impurity model Eq. (^) equals the total electron 
density {n) of the original lattice problem Eq. (0). 

< dUa > +Y.< Co,aCo,<x >= n (8) 

cr fj 

To determine Ud within the impurity model Eq. (P), we first note that the (i- level is 
renormalized from its bare value. This happens due to two sources, a short time Hartree-like 
renormalization due to the interactions and a long time contribution due to the hybridization 
which is described by the RG flow of section 0. 

To study the short time renormalization, we consider the non-fiipping part of the Hamil- 
tonian given in Appendix A. The shift of the d-level AEd comes from the change of the 
ground state energy of the conduction electron sea, when a (i-electron is present. The renor- 
malized (i-level is given by 

Edi^i) = e°^-fi + AEdifi) (9) 

Here AEd depends on the chemical potential /i since the shift of the ground state energy 
changes as /i is varied. 

There exists a critical chemical potential, fi = fi°, at which the renormalized d level lies 
at the chemical potential, 

EM = (10) 

In the absence of hybridization, such a condition indicates that, at /i = the heavy 
electron lies right at the Fermi level of the electron bath. Both Ud and Uc jump as 
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varies through /ic. Therefore, a finite range of electron densities correspond to the chemical 
potential fi°. Over this range of the electron densities, the d-level is pinned at the 
chemical potential. Within this pinned density range, there is no energy barrier for the 
transition from one local charge state to another. 

The essential question is whether, when the hybridization term is present, the chemical 
potential can still be adjusted to a value fic so that the renormalized d-level is equal to zero. 
In the language of critical phenomena, our procedure here is equivalent to renormalizing a 
massless field theory. In Section 0, we will show through the RG analysis that this can 
indeed be achieved when the hybridization is irrelevant. 



III. PARTITION FUNCTION OF THE GENERALIZED ANDERSON MODEL 

We now turn to the RG analysis of the generalized Anderson model Eq. (P) near the 
mixed valence regime where the renormalized d level is near the chemical potential. 

To derive the RG flow, our strategy is to first construct an expansion around the atomic 
limit, in terms of the hopping amplitudes between the local atomic configurations. It leads to 
a partition function written as a summation over histories associated with the local degrees 
of freedom. This in turn is equivalent to a 0+1 dimensional statistical mechanical model with 
long range interactions between the flipping events. Such an atomic expansion originates in 
the work of Anderson and Yuval in the Kondo problem and of Haldane in the asymmetric 
Anderson model [|1^. Here we generalize it to arbitrary local configurations, and construct 
a closed set of RG equations. We stress that, this atomic expansion is perturbative in the 
hoppings, but non-perturbative in the interactions and the symmetry breaking fields. 

To illustrate the methodology, we recall that, RG equations for Kondo-like impurity 



problems have been derived using a) poor man's scaling p0| , p!3[] ; b) multiphcative RG 
and c) mapping to one-dimensional statistical mechanical problem. |3y|T^ In the most gen- 
eral form, both the poor man's scaling and the multiplicative RG are difficult to carry out 
systematically: many couplings are present and it is not a priori clear which combinations 
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of these couplings come into the flow. We will show that, for a general set of local config- 
urations, the partition functions derived from the atomic expansion can be mapped onto a 
class of discrete classical spin chain models with a long range interaction in symmetry 
breaking fields. This will allow us to carry out a systematic RG analysis in the general case, 
by generalizing Cardy's analysis of these models PT| to incorporate the symmetry breaking 



fields. We note that, Haldane's work for the asymmetric Anderson model emphasizes 
the importance of a symmetry breaking field (i.e. the impurity level) as a consequence of 
the lack of particle-hole symmetry. Since the focus is on identifying various regimes within 
the strong coupling Fermi liquid phase, scaling equations are constructed for the hybridiza- 
tion and the symmetry breaking field only. In Cardy's work |^ for the general class of 
one dimensional classical discrete spin chains with long range interactions, RG flow is con- 
structed for all coupling constants, but without incorporating the symmetry breaking fields. 
The effect of the symmetry breaking fields in this scheme has been studied previously in 
other context. For our purposes of identifying different low energy fixed points and the 
associated transitions/crossovers between them, it is essential to study the renormalization 
of all the coupling constants in the generalized Anderson model. Meanwhile, particle-hole 
asymmetry generates non-zero effective fields even if we start from zero values of the bare 
fields. By combinging the formalism used in the works of Haldane and Cardy, we are able 
to derive the systematic RG equations for the generalized Anderson model near criticality, 
and identify and study Fermi liquid and non-Fermi liquid phases. 

In this section, we give the details of the atomic expansion for the partition function of 
the generalized Anderson model. In the next section, we derive the RG flow and establish 
the nature of the resulting fixed points. 

A. The Partition Function as a Summation over Histories 



For the generalized Anderson model at ?7 = cxd, the local degrees of freedom are 



characterized by the + 1 d-electron states, \a >= |0 >, or |(Jm > for m = 1, A^. 
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To perform the atomic expansion, we first separate the impurity Hamiltonian into two 
parts 



H = Ho + Hf (11) 
where Hq is diagonal in the space of the local states: 

k a c^a' 

+^ E ^VA'd.' + ^ E 4c.4^. (12) 

(T,<t' c 

The remaining part 

^/ = E^(4c. + /i-c.) + ^ E dld..c\,c, (13) 

flips from one local state to another. 

We can write Hq in terms of the projection operators X^a — \ol >< a\, 

H^^Y^HaX^a (14) 

a 

where 

H^^E^ + Y V2c\c^ + E Ekci^Ck^ (15) 

7 7 

Here cl^ creats a Wanniar state of the conduction electrons of spin 7 at the impurity site. 
The local levels are Eq — 0, Ea- — E^, and the potentials are given by 



— + — 

N N 

TV = (16) 



Meanwhile, the flipping part can be decomposed as follows. 



Hf=YQic^,f3) (17) 



where 
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Q{a,P) = \a><a\H0>< I3\ (18) 

describes the process flipping from the local state \13 > into \a >. Specifically, 

Q{(T,0)=Q\0,(T) = tdlc^ 
Y± 

Qia, a') = -j^did^.clc^{l - 5^,,) (19) 

The partition function 

Z = J DcDd exp[-{So + J^^ drHfir))] (20) 



can be expanded in Hf 

n=0 

where the transition amplitude, for a given initial (and final) state \a > and a sequence of 
flipping times ri < ... < r„, is given by 



r/S rn+1 _ 

^ = 2^ / dTn-: dTi... c/ri VA(a;r„,...,ri) (21) 
ZZ^ Jo Jo Jo ~7 



A(a;r„,...,ri) = (-1)"/ Dc Dd < a\T[expi-l3Ho)Hf{Tn)...Hf{n)...Hf{n)]\a > (22) 

For each transition amplitude, the path integral over the c?-states can be carried out 
through inserting a complete set of local states at every discrete imaginary time, leading to 

A{a; r„, n) = ^ Dc exp[-Hail3 - r„)]Q'(a, a„)... x 

02, ■■■."71 

xexp[-Ha,,^^{ri+i - rj)]Q'(ai+i, a^) exp[-Ha^{Ti - r^.i)]... x 
xexp[-Ha,,{T2 - n)] Q'{a2,a) exp[-HaTi] (23) 

Here, a, ^2, and ri,...,r„ label a Feynman trajectory: the local configuration starts 
from \ai >= \a >, changes from \ai > to > at time r^, (i = 1, ...,n — 1), and returns 
to >= I a > at r„. This is illustrated in Fig. 

The flipping operator (5'(«i+i)«i) =< > can be separated as follows, 

<a\Hj\P>=y'0\a,P) (24) 
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Here is the hopping amphtude associated with the flipping event (a,/3): ?/o,o- = Uafi—^^ 
y'^^, = -jj-. 0'{a,P) is composed of conduction electron operators which describe the dis- 
tortion in the sea of conduction electrons associated with the flipping event: 0'(cr, 0) = 
O'\0,a) = Co-, and 0'{a,a') = cl.,Co. It is the reaction of the conduction electron bath 
to these local distortions which will renormalize the hopping amplitude as we go to lower 
energies. 

We can now trace out the conduction electron degrees of freedom. The complication in 
evaluating each transition amplitude lies in the fact that, it has a history dependent potential. 
This can be dealt with most conveniently through the bosonization procedure. ||39|j3|] The 
relevant details of the bosonization applied to the impurity problems is summarized in 
Appendix A. For our impurity problem, we need to retain only the S-wave component of 
the conduction electrons. The conduction electron operator is represented by 

c.(x) = ^e-^*^(^) (25) 

where the $(j held is deflned in Appendix A. The radial dimension x has been extended 
to the negative half-axis. Hence, only one chiral component of the Tomonaga boson is 
introduced. 

The projected Hamiltonian transforms into, 

i^a=i^c + K + E— (^)-O (26) 

where po is the bare conduction electron density of states at the Fermi level. The phase 
shifts 62. are determined by the scattering potentials deflned in Eq. ([T6|), 

62 = 6i + 626^0 

S'o = (27) 



where 



= tan-^{7fpoVi/N) 
^2 = tan-\7ipoVl^/N) (28) 
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All the short time dynamics (for r smaller than ^o) are included in the shift of the atomic 
level, leading to the renormalized levels E'^ = + AE^,. The details for this level shift 
are given in Appendix A. The history dependent potential is treated through introducing a 
canonical transformation at each imaginary time 

Us = exp{i5^) (29) 

The transformed potential is time independent, due to the following property: 

UjH^Us = H, + —A^^o (30) 

7^ Po CiCC 



It is important to note that, the unitary transformation operator given in Eq. (p9D has 
the same form as the fermion operators. The transition amplitude now reduces to |^ 



a„+i=ai=a,a2,...,a„ 

n-1 

xexp[-E'^{n -rn)-Yl K^+ii'^i+i - n)] X 

i=2 



X < 0{an+i,an)iTn)...0{ai+i,ai){Ti)...0{a2,ai){Ti) > (31) 



Here 



0(ai+i, ai){Ti) = exp{HcTi)0{ai+i, ai)exp{-HcTi) (32) 

where the disorder operators are given by 

0(«,+i,«,) = iU^Usi^JO'ia,+i,a,){n^Ul) (33) 

It is the product of the canonical transformation operators, which represent the physics of 
a time dependent potential, and the conduction electron operators reflecting the distortion 
in the conduction electron sea. In terms of the boson fields, 

0(0, a) = exp{i{l - 62/TT)^a)exp{-iJ2i^i/'^)^j) 

7 

0(a, a') = exp(i(l - 62/n)^a')exp{-t{l - 62/71)^^) (34) 
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Eqs. ( p4D specify the initial values for the anomalous dimension e^^^ of the disorder 
operator 



Specifically, 



0{a,(3)^exp{tY^elp%) (35) 



TT 71 

ely = -(1 - + (1 - (36) 

TT TT 



Combining Eqs. ([21|) and (|3lD, we arrive at 

-E^:., / , "1 (37) 

j So 

where the transition amplitude is given by 

A{Tn,...,Ti) = < 0(a„+i,a„)(r„)...0(ai+i,ai)(ri)...0(a2,ai)(n) > 
= ^ / exp{-S,+ f Dr'Y.j^^{r')%{r')) 

lie J "'0 m 

= exv{- [ Dr'Dr"J2h{r')hir") < $ ,(r')<f ,(r") >) (38) 
with the "source current" 

hir)=T.:jir-r,)el^^^^^ (39) 

In Eq. (^7]), the cutoff is explicitly written in the integration range. Also the dimension- 
less form of the flipping amplitude has been introduced: Uat+ia, = u'ai+^aAo- Specifically, the 
charge and spin fugacities are given as follows, 

1/(0, (y) = yt = 

y{a ^ a') = y, = (40) 
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In deriving Eq. (pT)), we have also absorbed in Zq, besides the free conduction electron 
partition function Zc, a shift in the ground state energy such that J2a — 0- Specifically, 

N 



(iV + 1) 



Using the ^ behavior of the long time correlation function for the Tomonaga bosons, 
we finally arrive at 

where 



n- 



- E ^^(?/«.«... ) + E ^^^±^ (43) 

The partition function is now a summation over all possible histories of the local degrees 
of freedom which fluctuate between the + 1 local states \a >. Each history, labeled by 
{«!, an', Ti, Tn}, is a sequcncc of transitions between the local states from ai through 
a„ taking place at the time ri < ... < r„, as is illustrated in Fig. |1|. The action Eq. ( ^31) 
gives the statistical weight of such a history. 

We can interpret Eq. (^) as that associated with the partition function of a plasma of 
kinks with logarithmic interactions. It has (multi-component) charge e^^ and fugacity y^^. 



B. Equivalence to a Discrete Spin Chain with Long Range Interactions 

The partition function now has its action expressed as a summation over logarithmic 
interactions between the pair-wise kinks. Each kink is a flip in the original 'spin' space (i.e., 
the local states |a; >). The logarithmic interactions between the kink-pairs can therefore be 
transformed into interactions between the "spins". It is straightforward to show that, the 
action in Eq. ( ^3|) can be written in a symmetric form, 
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S[ti, ...Tn] = "i) + -^("j+i, "i+i) - K{ai, Oj+i) - K{ai+i, aj))ln 

So 



-T.Hya.a.,.)+T.ha j^'^y'^ (44) 

i i 

which can then be transformed into a long range interaction between the spins, 
S[ti, ...Tn] = K{ai, aj) ° 



Tj - Tif 



-T.Hya.a.,J + EK,. , (45) 

i i 

where the spins \a > can be in + 1 components. 

This transformed action is a special case of a class of discrete spin models with 
interactions considered by Cardy, ||2T| generalized to incorporate an effective field. The bare 



values of the stiffness constants —K{a,P) are determined by the bare values of the charges 



^(«,/^) = -^E«/3)' (46) 



2 7 



Specifically, 



KiO,a) = -et 

X(a,(T') = -e,(l-5.,.') (47) 



where the bare values of the charge and spin stiffness constants are 



2 TT n TT 

e° = (1 - ^)2 (48) 



IV. RG FLOW OF THE GENERALIZED ANDERSON MODEL 
AND THE NATURE OF THE FIXED POINTS 

In this section, we derive the RG equations using the atomic representation of the par- 
tition function derived in the previous section. The RG flow has the most general form 
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in terms of the charge and spin fugacities, the charge and spin stiffness constants, and the 
(i-level (a symmetry breaking field). The flow equations enable us to identify all the fixed 
points. They correspond to a universal strong coupling fixed point, a plane of weak coupling 
fixed points, and a plane of intermediate coupling fixed points. All these fixed points are 
studied in detail. We will establish in the next section that these fixed points describe Fermi 
liquid and non-Fermi liquid phases of the extended Hubbard model. 



A. Renormalization Group Equations 



We construct RG equations for small Ed in this section. As was discussed in Section |II C 
in certain range of parameters, the extended Hubbard model in infinite dimensions maps 
onto an impurity model which has a zero renormalized mass {d level E^). Therefore, the 
small Ed assumption is justified. However, since the RG procedure we use is perturbative in 
fugacities but non-perturbative in stiffness constants and fields, the RG equations for finite 
fields can also be constructed and are given in the Appendix B. 

The RG equations describe the flow of the dimensionless couplings as the bandwidth 
1/^ is reduced. The details of the derivations are presented in Appendix B, where arbitrary 
local states with a partition function of the form given in Eq. (^2| - ^5l) are considered. For 
our + 1 local state problem, the RG charges are the charge and spin fugacities i/t and yj 
(defined in Eq. (|40|) ), the charge and spin stiffness constants et and ej (defined in Eqs. ( ^Tf ) 
and (^8D ), and the dimensionless d level E^,^ = — h^. The RG equations are given as 
follows, 



dyt 



[I - tt)yt + {N- l)ytyj 



dln^ 
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dln^ 

n - -^yl - ^T/, (49) 



where F is the free energy. |H2 



In the renormahzation of the fugacities, the hnear terms give the associated anomalous 
dimensions, while the quadratic terms reflect the non-abelian nature of our + 1 state 
problem. Physically, these cross terms reflect the coupling between the local spin and charge 
degrees of freedom. The renormahzation of the stiffness constants reflect the correction to 
the interactions due to the fugacities of the charge and spin kinks. 

In the renormahzation of the energy level, the term arises because in integrating out 
the degrees of freedom the empty state gains a hybridization energy times as much as 
a spin state does. This reflects the particle-hole asymmetry of the mixed valence problem. 
131 On the other hand, a spin state gains exchange energy through the other (A^ — 1) spin 



states, leading to the term. We emphasize that, the symmetry breaking field is a relevant 
perturbation. 

For N=l, the above RG equations reduce to those of the anisotropic Kondo problem, 
through identifying t with with (1 — 2/-ntan^^{J\\^o/'^)Y, and —Ed with the magnetic 

field H applied on the local moment. [|1^,0 For N=2 and in the limit Vi = V2 = 0, we 



reproduce equations of Haldane's for the flow of the d level and the hybridization for 



the asymmetric Anderson model in the mixed valence regime. 

B. Fixed Points and Critical Behavior 

The kinks are bound together when the amount of attraction between them are strong 
enough. This leads to fugacities renormalized to zero. As in the Kosterlitz-Thouless transi- 



tion in the XY magnet, [^] unbinding of the kinks occur when the attractions between the 
kinks become weak. In our case, et and reflect the strength of the long range interactions 
between the kinks. Since the spin and charge kinks are coupled, the unbinding of these kinks 
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can occur either at the same time or at different stages. The above RG equations lead to 
three kinds of fixed points. The phase diagram in the — ej plane is given in Fig. |^. 

The weak coupling fixed points occur when both the spin and charge kinks are bound, 
corresponding to the hybridization and the exchange coupling both being irrelevant. This 
occurs when both the renormalized charge and spin stiffness constants are large enough, 

e*>l 

e* > 1 (50) 

where the stiffness constants are renormalized from the bare values by the fugacities. To 
second order in the fugacities, 

,*-,o_ , 2(iV + lK-(iV-l)e° _ , (iV-l)e° 
~ * 2(e?-l) ^y^' 2(e°-l) 

= - iy^rzA-, - A^(^J)^-^ (51) 

For N=l, this reduces to Si/n < -(V2 - 1) - t/T. PJig For AT > 2, it corresponds to 



iV(^f -2(l-^)^-(2-(l-^)2)>0 

^ < (52) 

TT 

The strong coupling mixed valence fixed point occurs when both the spin and charge 
kinks are unbound, corresponding to the hybridization and the exchange coupling both 
being relevant. The basin of attraction is given by 

el < 1 (53) 

and a range of e* < 1 when ej" > 1. 

This strong coupling fixed point is beyond the reach of perturbative RG. It is however in 
the same universality class as the usual strong coupling phase of the Anderson model: below 
a non-zero renormalized Fermi energy, both the local spin and charge degrees of freedom are 
quenched through the formation of quasiparticles. 
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In the spinful problem (A^ > 2), for a range of couplings within the following domain 



4>1 

e* < 1 (54) 

an intermediate phase can occur for which the spin kinks are unbound while the charge kinks 
are bound. Here the hybridization is irrelevant, while the exchange coupling is relevant. The 
local spins are quenched below a non-zero coherence energy, while the local charges are not. 
We note that, the RG equations imply that, when Uj becomes very large, it will start to 
drive yt to increase. However, this late stage is beyond the reach of perturbative RG. We 



will establish in Section |IVE| that, the fixed point with irrelevant yt and relevant yj is indeed 
stable. The transition to the strong coupling phase occurs when the effective charge stiffness 
constants, further renormalized due to the unbound spin-kink plasma, reaches the critical 
value 1. 

We emphasize that, from the quadratic couplings in the scaling equations, a relevant 
charge coupling necessarily drives the spin coupling to be relevant, even with a bare tj > 1. 
Physically, a spin kink (created by dj^d^) is the composite of two charge kinks (created by 
dcr). Unbounded charge kinks will strongly screen the interactions between the spin kinks, 
leading to the unbinding of the latter as well. 

The phase boundaries are given in Fig. ^ The vertical thick line separates the weak 
coupling phase from the strong coupling phase. The horizontal thick line separates the 
weak coupling phase from the intermediate coupling phase. The (schematic) dashed line 
represents the boundary between the intermediate coupling and strong coupling phases. 

We note the analogy of these three phases with those in the two dimensional defect- 



driven melting problem. |^ There, the binding-unbinding transitions are associated with 
the dislocation and the disclination defects. A dislocation is the composite of a pair of 
disclinations, in analogy to the spin kink being a composite of a pair of charge kinks. The 
hexatic phase exhibits bound disclinations and unbound dislocations, in analogy to our 
intermediate coupling phase with irrelevant charge fugacity and relevant spin fugacity. 
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In terms of the physical interactions, i.e. the density- density interaction Vi and the ex- 
change interaction V2, the basin of attraction of the weak couphng fixed points corresponds 
to a range of attractive density-density and ferromagnetic exchange interactions, while the 
intermediate coupled fixed points occur over a range of attractive density density and anti- 
ferromagnetic Kondo exchange interactions. We stress that, these are effective interactions 
at the beginning of the scaling trajectory. In generic models in infinite dimensions with a 
non-Lorentzian density of state, there are high energy dynamics which have to be integrated 
out before the low energy scaling regime is reached. We will discuss the implications of our 



results for more realistic models in Section VI 



We end this subsection with a discussion of the critical behavior of the zero temperature 
quantum phase transitions between these various states we identified. For reasons to be 
described at the end, these zero temperature quantum phase transitions can be formulated 
only in infinite dimensions. At finite dimensions we expect that the system will have sharp 
crossovers. 

The transitions are characterized by the collapse of an energy scale. 



1 



ej? ~ eoe (^--^c)'' (55) 

The transition between the weak coupling and strong coupling mixed valence states takes 
place when we vary the interactions (and hence ej' and e*) through the vertical thick line 
in Fig. |[ This one stage transition is illustrated in Fig. 0(a), where ei labels a line in the 
et — tj plane that passes through the vertical thick line at ec- The RG equations imply that, 
the unbinding transition for spin kinks is driven by that of the charge kinks. The critical 
behavior is therefore determined by the renormalization of yt only. Hence, ei;' in Eq. ( |55D 
corresponds to the Fermi energy of the Fermi liquid in the strong coupling state, and the 
exponent Tj = ^. 

When we vary the interactions (and hence et and ej) through the horizontal thick line 
and the dashed line in Fig. 0, two stages of transition take place, from the weak coupling 
through the intermediate coupling to the strong coupling phase. This is illustrated in Fig. 
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^(b), where €2 labels a line in the et — ej plane that passes through the horizontal thick 
line at ec2, and through the dashed line at ed- Near ec2, the ep in Eq. (|^) refers to the 
coherence energy for coherent spin excitations in the intermediate coupling state. Close to 
eel, on the other hand, ep refers to the coherence energy for the coherent charge excitations 
in the strong coupling phase. The exponent characterizing the unbinding of spin kinks at 
ec2 is given by 77 = 1. Due to rotational invariance, it differs from the exponents in the 
corresponding A^-component Coulomb gas problem, and is independent of N. The exponent 
characterizing the unbinding of charge kinks at ed, on the other hand, is r] = |. 

In the following subsections, we discuss in detail the nature of these individual phases. 



C. Strong Coupling Mixed Valence State 

Within the strong coupling mixed valence state, both the hybridization and the exchange 
coupling are relevant. The effect of interactions becomes smaller as we go to lower energies. 
The nature of the fixed point is that of the strong coupling mixed valence state of the usual 
Anderson model. The ground state is well described by various variational wavef unctions, 
such as those of Varma and Yafet. [^] The physics of this state can be systematically 



studied in the slave boson condensed phase in the large N approach. [^] Within the strong 



coupling phase, there exist the empty orbital, local moment and mixed valence regimes. 
||T3| , |3^ Crossover between these various regimes occur as the renormalized d-level is varied 
with respect to the Fermi level. 

In the following, we use our scaling equations to determine a generalized form for this 
d-level shift and resonance width A(^) = 7cyf{^)/^, in the presence of the finite density- 
density and exchange interactions. In the generic case, the flow of the d-\evel can only be 
derived through integrating numerically the closed set of the RG flow. For a vanishing bare 
spin fugacity, yj = 0, and when the renormalization of the stiffness constants is ignored, the 
leading terms for the running d-Wel and the running resonance width can be determined 
from 
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J^Ml-2e,)A4)-.-" (56) 



where Aq = n{y°Y/$,o- This leads to 

A 1 - )-(2'^'-i) 



TT 2ei - 1 

A(0=Ao(e/eo)-('^'-') (57) 



from which we construct the invariant level and resonance width 



vr 2et-l 

A*= Ao(Ao^o)^^ (58) 
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In terms of these scaling invariants, the running d level and resonance width are given by 

A* 1 - ('£A*)"(2et-l) 

m)= e: + {n- 1) — 

A(0= A*(eA*)-(2^*-^) (59) 

A logarithmic shift is recovered when = = 1/2 is used, corresponding to Vi = V2 = 0. 

At low temperatures, the crossover between the various regimes can be determined from 
comparing the invariant d level and the invariant resonance width A*. The local moment 
and empty orbital regimes occur when << —A* and E^ » A* respectively. For 
1^*1 ^ ^Yie system is in the mixed valence regime. This crossover behavior is illustrated 
in Fig. §(a). Here we have used Eq. (^, so that the variation of E^ is described in terms 
of the chemical potential fi of the lattice model. At zero temperature, the mixed valence 
crossover extends over a scale of ~ A*. This crossover can also be illustrated through the n 
vs fi curve, shown schematically in Fig. ^(a) and to be discussed further in Section [V|. 

Finite temperature (or frequency) cuts off the scaling. As temperature is increased, both 
the running (i- level Ed{T) and the running resonance width A(T) vary. The crossover tem- 
perature scale, separating the high temperature mixed valence regime from the low tempera- 
ture empty orbital or local moment regime, can be determined from Ed{T) = max(T, A(T)). 
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This specifies the crossover hnes— the dashed hnes— in Fig. §(a). In the local moment 
regime, there will be a further crossover at the scale of the Kondo temperature (not shown 
in the figure). In all these regimes, the low energy behavior is described by the Fermi liquid 
theory. The renormalized Fermi energy is given by the renormalized resonance width A*. 



D. Weak Coupling Mixed Valence State 

At the weak coupling fixed points, the fugacities become smaller as we go down to 
lower energies. Therefore, the atomic expansion in terms of the running fugacities become 
more justified at lower energies. All the correlation functions can be derived from the RG 
procedure. 

The running fugacities and d level can be derived from integrating the RG flow equations: 

E,{oo) = El- (60) 

So 

where A = ^'^~P^J^ , and E*, = E'. + ^. In this case, the level shift due to the particle- hole 
asymmetry is small, since the hybridization is irrelevant. 

The crossover diagram is given in Fig. §(b). There again exist local moment, empty 
orbital, and mixed valence regimes. The important difference with the the crossover diagram 
of the strong coupling phase (Fig. Ka)) is that, in this case the renormalized resonance width 
A* vanishes at zero temperature; A zero temperature critical point occurs at yU = yUc- This 
can also be shown through the n vs fi curve given in Fig. |^(b), where n is the total electron 
density. 

To determine the finite temperature crossovers between the different regimes, we define 
the running "resonance width" A(T) defined by 

A(r) = vry^T = nt%{TQ^'<-'^ (61) 
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The local moment, empty orbital, and the mixed valence regimes occur when Ed{T) « 
-max(T, A(T)), Ed{T) » max(T,A(T)), and |^d(OI < max(A(T) respectively. This 
determines the crossover behavior illustrated in Figs, ^(b) and ^(b). 

For ~ (— A/^o) < <~ v4/^o, an additional crossover occurs at an energy scale T', 

T' = -(^^)^ (62) 

which is determined from \Ed{T')\ ~ A(T'). The system is in the critical regime at temper- 
atures higher than T', since here it is essentially an Ed{T) = problem, with a temperatur- 
dependent resonance width which is smaller than the temperature. In particular, at E^ = 0, 
i.e. /i = {ic, this mixed valence critical regime persists all the way down to zero temperature. 
The mixed valence state corresponds to a massless critical point. 

The mixed valence phase here is characteristically different from the strong coupling 
mixed valence phase. The irrelevant hybridization makes the running resonance width to be 
always smaller than the scale, A(u;) < u. Local correlation functions in the mixed valence 



regime can be calculated through the RG procedure. p2 

Here we focus on the local Green's functions, which is a 2 x 2 matrix for each spin 
component: 



G 



( n n ^ 



y Gdc Gdd J 



(63) 



The details of the RG calculation for the renormalization of the exponents can be found 
in Appendix C. The results are as follows. Gdc and Gdd exhibit infrared divergences, and 
have the following low energy behavior, 

Gdcico) ~ 

Gddiu;) ~ ^-1+^ (64) 
Here the exponents are, to leading order, as follows 

a = -{61 + 6;)/n + {{61 + 6l)/7,f + {N - l){5l/ixf 
(3={{6l + 6l)l7:f + {N-l){6ll7^f (65) 
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where 5*2 are the renormahzed phase shifts determined from Eq. (0). The algebraic 
divergences in Gdc and Gdd reflect the fact that, both the impurity spin and charge degrees 
of the freedom are not quenched. 

On the other hand, the local c-electron Green's function Gcc does not show infrared 
divergences. Formally, creating a local c-electron does not creat any kink in the path integral 
representation. Hence, Gcc is not renormahzed from its fixed point form. At the mixed 
valence fixed point, Gcc has a long time ^ behavior. The absence of infrared divergences 
in the electron bath spectral function will be used in Sec. [VD| to argue that, our results 
derived using a Lorentzian density of states can be applied to the case of generic density of 
states as well. 



E. Intermediate Coupling Mixed Valence State 

In this case, the RG equations yield an initially decreasing charge fugacity yt and an 
initially increasing spin fugacity Uj. When yj grows to order unity, the flow equations would 
indicate that yt starts to grow. However, this is beyond the regime of validity for the 
perturbative RG. 

To understand qualitatively the nature of the fixed point, we analyze the model in the 
limit that the antiferromagnetic exchange V2 is much larger than the conduction electron 
band width, while the hybridization t is much smaller than the band width. This serves as 
an effective Hamiltonian that the original Hamiltonian renormalizes to at an intermediate 
scale. 

U ^ 

Heff = ^ did^>cl^,Co,a + ViY, dld„cl^^,Co,a' + Crf I] dld„ + 1^Y.Y. 4„diadl^,dia' 

a,a' (t,(t' f i u^a' 

+ ^t(4co,a + /i-C.) + Xl^l(Cl,<7Co,a + ^-C.) + Y tij{cl^Cj^a + h.C.) (66) 

Here c| ^ creats effective neighboring states which are coupled to the conduction electron 
state at the impurity site that cj creats. We consider N=2 for simplicity. The parameters 
are such that V2 » Vi,ti » t, and U = 00. 
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To order 0(|^■,■^|■), diagonalizing the V2 coupling first leads to singlet and triplet sub- 
sectors within the = 1 sector. And the triplet sector locates at high energies. The mixed 
valence condition is satisfied at an e'^ such that the singlet sector and the = sector are 
(nearly) degenerate. We can now integrate out the high energy states. This leads to the 
following effective Hamiltonian for the low energy sectors, 

k 

+ J2 (''^'k,k's'^eck,^ck,^asgn{a) + h.c.) + J2 ^'A;,fc'(s"^s - e^e)4,<xCfc,a (67) 

k,k',cr k,k',a 

Here and creat the local singlet and the local empty state respectively, while 4 ^ creats 
conduction electron scattering states. The effective interaction and the effective "hybridiza- 
tion" amplitude are given as follows 



t'k,k'-^ (68) 

Given that t is small, and < 0, the effective "hybridization" is marginally irrelevant. 
Therefore, the spin and charge sectors are indeed decoupled asymptotically. 

The crossover diagram is also given qualitatively by Fig. Kb), except for the local 
moment regime in which the spin degrees of freedom is eventually quenched through the 
Kondo process. The n vs. /i curve at zero temperature is again qualitatively given by Fig. 
1(b). 

For the mixed valence regime, the spectral functions have the form of a convolution 
between the resonance in the spin sector and the algebraic term in the charge sector, and 
therefore also exhibit an algebraic behavior at low energies. The exponents can be deter- 
mined from the effective charge stiffness constants, renormalized (further from Cj) due to 
the screening of the free spin kink plasma. This further renormalization can be analyzed 
in detail by following the scaling trajectory and using Eq. (|67D, in the same spirit as the 



analysis of the hexatic phase using the Debye-Hiickel approximation. ||3^ In particular, the 
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exponent is universal when the critical line separating the intermediate coupling and the 
strong coupling states, i.e. the dashed line in Fig. 0, is approached. 



V. NON-FERMI LIQUID AND FERMI LIQUID PHASES 
OF THE EXTENDED HUBBARD MODEL 

In this section, we establish the implications of the results of the associated impurity prob- 
lem in the previous two sections on the phases in the extended Hubbard model. We show 
that, the phases corresponding to the weak coupling and intermediate coupling mixed va- 
lence states are metallic non- Fermi liquid phases. They have vanishing quasiparticle residue, 
algebraic local correlation functions, and asymptotically separated charge and spin excita- 
tions. The phase corresponding to the strong coupling mixed valence state is a usual Fermi 
liquid. 



A. Metallic non-Fermi Liquid: Weak Coupling 
and Intermediate Coupling Phases 

Consider first the phase corresponding to the weak coupling mixed valence state. We first 
establish that, the critical chemical potential fic, at which the mixed valence state persists 
to zero temperature, corresponds to a range of electron densities. 

The local correlation functions of the extended Hubbard model in infinite dimensions 
are given by the impurity problem. In particular, the occupation numbers in the lattice, 
for a given chemical potential, can be obtained from the local Green's functions of the 
corresponding impurity model. Therefore, our analyses in Sections [II C| and [IV D| establish 



that, at zero temperature and also n = n^ + nc are discontinuous functions of the chemical 
potential: for fi > fic, ^ 1 — 0(t^) while for fi < fic, ^ O(t^). At finite temperatures, 
the discontinuity is smoothed out and changes into a fast crossover. Here we take fi = fic 
which again corresponds to a finite range of densities in the lattice model. 
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It is remarkable that the fact that our impurity model is associated to a lattice problem 
forces the impurity model to be exactly at criticality, with a larger symmetry than we would 
have naively expected. The essential point is that, in order for an incoherent state to be 
metaUic, it is necessary to allow for charge transfer between the localized degrees of freedom 
and the bath. This can only happen if the local charge degree of freedom is in equilibrium 
with the conduction electron bath. This requires the heavy level to be at the chemical 
potential. 

To determine the excitation spectrum of the system, we consider the lattice Green's func- 
tion. In the limit of infinite dimensions, all the lattice correlation functions can be recon- 
structed from the local correlation functions calculated in the associated impurity problems. 
The lattice Green's function. 



G(k,T) 



( 



\ 



Gai%T) 

GddiK t) j 

- < Tck(T)4(0) > 

- < T4(r)4(0) > 

can be determined from the Dyson equation for the lattice model 
/ 



V 



- < Tck(T)4(0) > 

- < T4(r)4(0) > 



(69) 



I 



(70) 
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where the momentum-independent self-energy can be calculated directly from the local 
Green's functions of the impurity model 
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This leads to, for A'" > 1, the self-energies with the following form at low energies. 
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1 7c. 



~ (o;)^^^ (72) 

Here the exponents jcc = 2e* — 2, 7(i(i = 1 — a, and = 2e* + P — 1, where a and /3 are the 
exponents associated with Gdd and Gdc given in Eq. (|65|). 

The form of the self-energies given in Eq. (|72[) destroys the quasiparticle pole. Therefore, 
the quasiparticle residue vanishes; all the single particle excitations are incoherent. 

We now turn to various physical correlation functions. Since both the impurity spin and 
the charge couplings are irrelevant, we expect that the multiparticle correlation functions in 
both impurity spin and charge channels show power law behavior. The spin and charge cor- 
relation functions are generally renormalized differently. The precise values of the exponents 
can be derived using RG as discussed in Appendix C. For example, 

< did,iT)di,d,iO) >~ (r)-"^ 
<4c^,(r)4,f/.(0) >~ (r)-"^ 

<44,(r)c^,t/.(0)>~(r)-°3 (73) 

where to leading order 

ai = 2e*(l-5^,,0 

a2 = 2e,X,,, +2e:^.(l-5,,,0 

as = 2e:(-^, -^)5.,.. + 2e:,.(-^, -^)(1 - <5.,.0 (74) 

71 71 71 71 

in which e*,.(f , f ) = ((1 - 5)2 + + f )2 + (iV - 2)(5)2)/2. 

These exponents correspond to a divergent rf^c^ pairing susceptibility and vanishing ex- 
citonic and d-electron correlation functions at low energies. The divergent pairing suscep- 
tibility signals a superconducting instability at a finite temperature Tc, with an anomalous 
normal state described by the non-Fermi liquid state. We will comment on this further in 



section VI 
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We now turn to the intermediate coupling phase. A critical chemical potential /ic also 
exists due to the irrelevant hybridization, as was established in Section [IV li| . It again 
corresponds to a finite range of densities in the lattice model, at which the mixed valence 
behavior persists to zero temperature. 

Since the exchange coupling is relevant, a resonance forms in the spin channel below 
the corresponding coherence energy. Therefore, spin excitations form quasiparticles which 
are expected to have Fermi-surface features. The local spin susceptibility is regular. On 
the other hand, the local charge susceptibility retains the algebraic behavior. The excitonic 
correlation function still vanishes in the zero frequency limit, while the superconducting Sc^ 
susceptibility still diverges at low frequencies algebraically: the system is again expected to 
be a superconductor with an anomalous normal state. 

To summarize, the non- Fermi liquid phases characterize incoherent metals with vanishing 
quasiparticle residue, self-similar local correlation functions, and asymptotically decoupled 
charge and spin excitations. In addition, they occur over a range of electron densities due 
to the pinning effect. 

B. Fermi Liquid Phase 

We now turn to the region where the hybridization is relevant. The system then flows 
away from the zero fugacity fixed points. The new fixed point is beyond the reach of the 
perturbative RG. The nature of the fixed point can be studied by other approaches describing 
the strong coupling phases. One approach is the large N limit of the model with fixed finite 
values of Vi and V2. 

In the large N mean field theory applied directly to the lattice model (|l|), the renormal- 
ized hybridization never vanishes, and the density is a continuous function of the chemical 
potential. The (i-electron spectral function is smeared into a coherent Abrikosov-Suhl reso- 
nance. Low energy excitations are coherent, and are composed of conduction electrons and 
heavy electrons with finite hybridization. The system is a Fermi-liquid. 
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That the large N mean field theory forces the system to be a Fermi liquid can be easily 
understood from our scaling equations (^). When the large N limit is taken, with a fixed 
finite values of Vi and V2, the associated phase shifts given in Eqs. ( pTf ) and ( pH]) vanish. 
Hence, the anomalous dimensions of the spin and charge fiippings, defined in Eqs. ( ^Hf ) 
and (0), fall in the basin of attraction of the strong coupling fixed points. Physically, the 
excitonic and orthogonality contributions are suppressed by a factor of 

C. Qualitative Phase Diagram of the Extended Hubbard model 

The qualitative phase diagram, in terms of the interaction strength and electron density, 
is summarized in Fig. ^ 

Fig. |^(a) corresponds to varying the interaction strength (in the — ej plane in Fig. 
0), as that characterized in Fig. ^(a). For interaction strength such that ei is above the 
threshold value ec given in Fig. 0(a), and over a range of electron density given by the solid- 
line shaded region, a mixed valence metallic non-Fermi liquid state occurs down to lowest 
temperatures (see, however, the discussion about long range ordering in Section [VT]). This 
range of density corresponds to the critical chemical potential fic in Figs. |l(b) and §(b). For 
electron densities outside this range, within the dashed-line shaded region, non-Fermi liquid 
behavior occurs above a crossover temperature, shown in Fig. §(b). The non- Fermi liquid 
state is characterized by incoherent charge and spin excitations, which are asymptotically 
decoupled. 

For interaction strength such that ei is below the threshold value ec, the system is a 
Fermi liquid. It corresponds to the strong coupling mixed valence state of the associated 
generalized Anderson model. Near the threshold interaction strength, the Fermi energy is 
given by Eq. ( |55D with the exponent 1] = ^. 

In Fig. |^(b), we vary the interaction strength as that characterized in Fig. ^(b). For 
interaction strength such that €2 is above the threshold value ec2 in Fig. 0(b), a (weak cou- 
pling mixed valence) metallic non-Fermi liquid phase again occurs down to zero temperature 
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over a range of electron density within the sohd-horizontal-hne shaded region, and above a 
crossover temperature for electron densities within the dashed-horizontal-line shaded region. 
The non- Fermi liquid state is the same as that in Fig. ^(a). 

For interaction strength such that €2 is between ed and ec2 of Fig. 0(b), and over 
a range of electron density given by the solid-vertical-line shaded region, an (intermediate 
coupling mixed valence) metallic non- Fermi liquid phase occurs down to lowest temperatures. 
For electron densities within the dashed-vertical-line shaded region, the non-Fermi liquid 
behavior occurs above a crossover temperature given in Fig. H(b). The non- Fermi liquid 
state is characterized by coherent spin excitations and incoherent charge excitations, which 
also decouple asymptotically at low energies. The spin coherence energy near is given by 
Eq. (|55D with the exponent 77 = 1. 

Finally, for interaction strength such that €2 is below the threshold value ed, the system 
is a Fermi liquid, corresponding to the strong coupling state. 



D. Dependence on the Form of the Density of States 

Thus far, we have focused on the extended Hubbard model with a Lorentzian density 
of states. The non-Fermi liquid phases occur due to the persistence of asymptotically de- 
generate states at the Fermi level. This occurs because the quantum hopping amplitude 
(the hybridization) which mixes these degenerate states is irrelevant. These couplings are 
irrelevant because of the slow relaxation of the electron bath to the hopping between differ- 
ent local charge states. Such a large reaction exists whenever there are non-zero density of 
states of the electron bath at the Fermi level. In the Lorentzian case, there is always a finite 
density of states for the bath electron. 

To establish the generality of our results, we consider the extended Hubbard model 
defined on the Bethe lattice with infinite connectivity. In this case, the bare density of 
states is bounded. P5|,p!0| The Weiss field G"^, which is a 2 x 2 matrix in our case, can be 



explicitly related to the full Green's function. 
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Go ^{iuJn)=go ^(iuJn) -^toiGijtjo (75) 

ij 

where go~^ is the inverse propogator of the corresponding non-interacting problem. More 
exphcitly, 

{G-'),, = iUn + fi-{tcc?G,, (76) 

Such a relation demonstrates that, when the renormalized c-electron density of states is 
non zero at the chemical potential, the Weiss field (G^ ^)cc will have a non-zero imaginary 
part. When the hybridization t is small, the c-electron acts as the electron bath that the 
local c?— electron degrees of freedom couple to. Sine the c— electron Green's function does not 



show infrared divergences, as is discussed in Section [IV D| and C, the local Green's function 
of the electron bath has the long time asymptotic ^ behavior. This establishes that, a finite 
density of states for the electron bath at the Fermi level is indeed self-consistent. This finite 
density of states can be effectively modeled as the inverse Lorentzian width. p6| 

Therefore, our results which are explicitly derived for the extended Hubbard model with 
a Lorentzian density of states is expected to reflect qualitatively the physics of the extended 
Hubbard model away from half-filling. At the quantitative level, one has to determine the 
charge and spin stiffness constants and e^, used in constructing the phase diagram in Fig. 
0, through integrating out non-universal high energy dynamics. The phase diagram will be 
altered correspondingly. 

Finally, within the model with zero direct ci— electron hopping, only (G~^)cc has a non- 
trivial form. Hence, direct (i— electron hopping effects are not generated. How new energy 
scales might be generated from an explicit direct (i— electron hopping term is an important 
issue to be further studied. 
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VI. CONCLUSION 



We have carried out a detailed study of an extended Hubbard model in infinite di- 
mensions, where the model reduces to a generalized asymmetric Anderson problem in a 
self-consistent medium. Within the associated impurity problem, we carried out a system- 
atic renormalization group analysis, and identified three kinds of mixed valence behavior. 
Each of these states is associated with a distinctive flow of spin and charge couplings. The 
usual strong coupling mixed valence state has spin and charge couplings both relevant; The 
weak coupling mixed valence state has spin and charge couplings both irrelevant; And the 
intermediate coupling mixed valence state has a relevant spin coupling and irrelevant charge 
coupling. The strong coupling state corresponds to a Fermi liquid phase of the extended 
Hubbard model, while both the weak coupling and the intermediate coupling mixed valence 
states describe metallic non-Fermi liquid phases of the extended Hubbard model. 

The non-Fermi liquid phases are metallic states with incoherent charge excitations. 
Within these incoherent metallic states, the local charge modes are necessarily close to 
the Fermi level for charge transport to occur. Such a pinning effect makes the non-Fermi 
liquid phases to occur over a range of electron densities. Hence, no fine tuning of parame- 
ters is necessary. The pinning of the chemical potential leads to an infinite compressibility 
at zero temperature: there is a jump in the n vs. /i curve, as shown in Fig. ^(b). It is 
interesting to note that, an infinite compressibility has been demonstrated in the numerical 



studies of the Hubbard model in the limit of small dopings away from half-filling. 147 



The non-Fermi liquid phases have vanishing quasiparticle residue: all excitations are 
incoherent. In addition, charge and spin excitations decouple asymptotically at low energies. 
Finally, certain local correlation functions show algebraic behavior, reflecting the self-similar 
nature of these incoherent metallic states. These are the local correlations in both the charge 
and spin channels in the weak coupling phase, and in charge but not spin channel in the 
intermediate coupling phase. These properties have striking analogies to the Luttinger liquid 
in one dimension. P,|3| 
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The novel phase that we identify, i.e. the intermediate couphng phase, has coherent spin 
excitations and incoherent charge excitations. In connection with the possible realization 
of such a phase in strongly correlated electron systems, we note that the normal state 
of the high Tc copper oxides exhibits low energy spin dynamics which are well described 
by renormalized quasiparticles with a Luttinger Fermi surface, but charge dynamics which 
appear to be more anomalous. [^] We will discuss physical properties of our intermediate 



phase in more detail elsewhere. 

We stress that at the end one ought to be able to express our results, derived using a 
local method, in a language that would make contact with perturbative calculations in the 
interactions. When this is done we may be able to identify some similarities between our 
results and those derived from other methods. This will also allow us to extend our results 
to models at finite dimensions. In this regard, the Fermi liquid phase is described in terms 
of the boson "condensed" phase within the slave boson formalism: the hybridization renor- 
malization is determined by the the quasi-condensate of the slave boson, and the ci-level 
renormalization is given by the mean field value of a lagrangian multiplier enforcing the no- 
double-occupancy constraint. [^] Our intermediate coupling incoherent metallic state has 
similarities with the boson non-condensed phase, which also has coherent spin excita- 
tion ("spinon") and incoherent charge excitation ("holon"). From our local treatment, the 
characteristics of an incoherent metallic state lies in its self-similar nature, a problem with- 
out scale. In both treatments of these incoherent metallic states, the local gauge invariance 
which is broken by the coherent hopping is finally restored. 

We now comment on several aspects of our results which are important to address in 
order to apply our results to finite dimensions as well as to more general models. 

In general, the limit of infinite dimensions separates the single particle and collective 
excitations. As a result, when the system undergoes an ordering instability, the one parti- 
cle property does not show precursor effect; The feedback of the long range order over the 
single particle properties is down by powers of ^. Within the incoherent metallic phases, 
there exist asymptotically degenerate local configurations, which is manifested in the alge- 
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braic low energy behavior in various local correlation functions. In finite dimensions, the 
coupling of the long range order on the local physics modifies the local effective action at 
temperatures below the onset temperature Tc, and has the effect of cutting off the power 
law singularities. In our incoherent metallic states, the most divergent local susceptibility is 
a pairing susceptibility. Therefore, the instability is expected to be in the superconducting 
channel. Our incoherent metallic states are then expected to be relevant at temperatures 
Tc << T << ujc in the same sense that the paramagnetic Mott insulating phase is well 
defined for T/v << T << Ci? (T/v being the Neel temperature). 

Within the extended Hubbard model in infinite dimensions with a Lorentzian density of 
states, the weak coupling incoherent metallic state occurs over a range of attractive density- 
density interactions and ferromagnetic spin interactions, while the intermediate coupling 
incoherent metallic state occurs over a range of attractive charge interactions and antiferro- 
magnetic spin interactions. We emphasize that, a Lorentzian form for the density of states 
suppresses the high energy dynamics, enforcing the system to be always in the scaling regime. 
In the generic case, the electron bath will be in the universal regime only within a narrow 
energy range around the Fermi level. The high energy dynamics have to be integrated out 
before the renormalization group results can apply. Our results here, derived for an extended 
Hubbard model with a Lorentzian density of states, give the classification of the universality 
classes for a general + 1 local configuration problem. To determine the basin of attrac- 
tion of the various phases, one has to first integrate out the high energy dynamics beyond 
the scaling regime. Such a procedure makes it possible that realistic Hamiltonians with all 
repulsive atomic interactions fall in the basin of attraction of the non-Fermi liquid phases. 
||50| Within generic models, the basin of attraction of non-trivial phases can be determined 
numerically. 

We conclude with a note on the methodology implications of our analysis. The general 
renormalization group scheme developed here opens a local view for strongly correlated 
electron systems. The universality classes of a correlated electron system can be determined 
by a small number of atomic configurations coupled to an electron bath, which can be treated 
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through the general atomic expansion plus the renormalization group procedure we outlined 
here. The specifications of the phase diagram, on the other hand, involves integrating out 
non-universal high energy dynamics, which falls in the domain of quantum chemistry. [pO|] 
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APPENDIX A: BOSONIZATION AND THE LEVEL SHIFT 



In this appendix, we summarize the bosonization procedure |3^J^ relevant to our discus- 
sion. We also discuss the shift in the atomic levels that should be incorporated in carrying 
out the bosonization. 

Consider the projected Hamiltonian given in Eq. (|T3|), 

H^ = E^ + Y. KcU^ + (Al) 

where 

H, = Y.Ekclck^ (A2) 

k,(7 

For the asymptotic behavior, as long as the density of states of the conduction electron is 
finite, a bosonization scheme can be used. For the purpose of studying an impurity problem, 
we need to consider only the S-wave component of the conduction electrons. We label the 
radial dimension by x, which is extended to the negative hafl-axis. 

The Tomonaga bosons are defined as follows. 



'27r 
L 



V 9-^ k 



bq,a = -ij ^ 4,aCfc+g,a (A3) 
^ 1^ k 

where q > 0, and L is the size of the radial dimension. 

The free electron Hamiltonian has the following form in terms of the Tomonaga bosons, 

H, = Y,VFqhlX,<^ (A4) 

q,<j 

where the boson dispersion is determined by vp = ^1 Po, the inverse of the density of states 
at the Fermi level. 

The conduction electron operator has the following form, 

c.(x) = -^e-^*-(^) (A5) 
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where 

^Ax) = E + b,e'n (A6) 

q \ 1-^ 

The electron density, Po-(x) = c^Cg-, is related to the $0- field as follows, 

P.{x) = (A7) 

Therefore, the potential term can be represented as follows, 

V:ctc^ = ^(^).=o (A8) 
irpo dx 

We use the prescription that, all quantities in the bosonization representation corresponds 
to the normal ordered quantities. In so doing, a Hartree-like shift in energy has to be 
incorporated separately. The exact form of such a level shift was discussed in detail for the 
spinless case in Ref. [^. It can be easily generalized to arbitrary spin degeneracy N, leading 
to, 

= ^ E EiHgo') - H9o' - V:)}e^-"'' (A9) 

where goiiujn) = Z^fc -^^ +/j-efc • "^^^ linear in V term in the above shift AE'q, is given by the 
g = component of the Tomonaga boson. 

Such a shift can be captured by expressing Ha in the bosonization form, 

Ha = H, + EJ + E ^(^).=o (AlO) 
^ vrpo dx 

where E'=Ea + AE^. 



APPENDIX B: RENORMALIZATION GROUP EQUATIONS 

In this appendix, we present the details of the derivation for the RG equations, which 
will be written in the general form, for M arbitrary local configurations \a >. 

The RG equations describe the fiow of the dimensionless couplings as the bandwidth 1/^ 
is reduced. The RG charges are the fugacities ya/s, the stiffness constants —K{a,P), and 
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the symmetry breaking fields ha- They satisfy yaa = 0, K{a, a) = 0, and J2a ha = 0. These 
relations are preserved in the renormalization process. 

Integrating out the degrees of freedom in the range {^,C, + dC,) leads to three kinds of 
renormalization contributions. One contribution arises from rescaling the old cutoff ^ by the 
new one ^' = ^ + in Eqs. (H2-H3p. This amounts to a renormalization of the fugacities 



and the fields, 

^ KXl + dlni) (Bl) 

The other contributions arise from the change in the integration range in the partition 
function, Eq. (^). The dominant renormalization arises from integrating out kink-pairs 
whose separation falls in the range {S,-,^, + dE,). There are two kinds of contributions in this 
category. 

The contributions from integrating out non-neutral pairs, i.e. pairs with the final state 
of the second kink different from the initial state of the first kink, amounts to a change in 
the fugacities: 

7 

for ai+i 7^ at. 

On the other hand, integrating out neutral pairs leads to a renormalization of the stiffness 
constants, 

7 

-dilnO E yl,,e'-^-''iKg^^„g^ + K^^,, - K^^^,,,) (B3) 

7 

as well as a renormalization in the field and in the free energy, 

7 7li72 



M 

71:72 
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Here the separation into the field renormahzation and free energy renormahzation is unam- 
biguously specified by requiring that J2a ha — he preserved in the process of renormahza- 
tion. 

Collecting all these terms, we derive the following general form of the RG equations. 

^^j^ = - E Vl.e'"-'^ (5) + K{a, 7) - K{(5, 7)) 
7 



APPENDIX C: RENORMALIZATION OF CORRELATION FUNCTIONS 
IN THE WEAK COUPLING REGIME 

In this appendix, we give the details of the derivation for the renormahzation of the 
exponents of the correlation functions in the weak coupling regime. 
The Green's function 

G^^{t) = - < TV'(t)V^(0) > (CI) 

where ■^^ = (d^, c^), can be expanded in terms of the hipping part of the Hamiltonian Hf, 

G^^r)^ ~ J Dc Dd T['ilj{r)ilj\0)exp[-{So + J^rHfirm 

^^E / dTn... dn... (iri /l(a;r;0;T„, ...,Ti) (C2) 

^ „=0 •'^o •'^o 

where the transition amplitude 

A{a; T- 0; r„, n) = (-1)" j DcDd 

< a\T[exp{-pHo)i;{r)Hf{rn)...Hf{n)...Hf{n)i;\0)]\a > (C3) 
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Without loss of generality, we take r > 0. For a given history similar to that illustrated 
in Fig. 0, the external time r lie between two flips, say t^-i < t < t^. The transition 
amplitude has the following form once a complete set of local states is inserted at every 
discrete time 



A{a] t; 0; r„, n) = (-1)" ^ Dc exp[-Ho,{(3 - r„)]Q(a, x 
xexp[-Ha^^^{Tk+i - Tk)] Q{ak+i, am)exp[-Ha^{Tk - r)] x 
xtp exp[~Ha^,{T - rk-i)]Q{am',ak-i) exp[-ifa^_^ (rfc_i - rfc_2)]--- x 
xexp[-Ha2{T2-Ti)] Q{a2,ai) exp[-HaiTi] (C4) 

where \am > and la^' > label states after and before the external time r. The operator ip'^ , 
on the other hand, flips from the local state \a > to \ai >. We emphasize that, the external 
fermion operators place constraints on the hopping history of the impurity: the insertion of 
a (i— electron creats a charge kink, while the insertion of a local c— electron does not creat 
any kink. 

Tracing out the conduction electron degrees of freedom, 
A{a; r; 0; r„, n) = Z, ^ ya„+u»„-y'a,+„aya^,,a,.,-y'a,,a, X 

Qn+i=Q,ai,...,a„ 
n-1 

XeXp[E'^Tn - E'^Jl - E'^^^^{Ti+i - X;)] < 0{an+l,an){Tn)--0{ak+l,am){Tk) X 

■t=2 

x*a„,a„,(r)O(a^/,afe_i)(rfc_i)...O(a2,ai)(ri)^t^^^(0) > (C5) 
where ^'^a^p = i^a,i3^ ^l,/?) have the following form, 

= exp{iel^^{T))6a,\o>Sp,\a> 

= exp{iel%{T))6o,,p (C6) 



with, 



el = -5^,, (C7) 
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The transition amplitude now has again a simple form in boson representation through 

< 0(Q;„+i,Q;„)(T„)...0(Q;jfc+i,Q;„)(Tjfc)*^(T)0(Q;^/,Q;jk_i)(Tjt_i)...0 > 

= J exp{-S, + Dt'Y.j,{t')%{t')) (C8) 



where 



=7 

''Oik,am 



- ri)el^,a, + S{r' - r„)e2„,„ + 6{t' - r)el + S{r')el, (C9) 

Here a primed sum means a summation which excludes n,l,k — l,k. Therefore, 

< 0(q;„+i, an){Tn)...0{ak+u oim){rk)^a{r)0{am': ak-i){Tk-i)...0{a2, q;i)(ti)*^(0) > 

^expi-^A'] (CIO) 

7 

where 

A J "'^ 7 7 1 '^hl I l'^^"'^^" y y ] a- ^^^^ ''"'^ ^ 

ii 12 '^'^ ii 12 

+E'E"eZ,,+i,a,,e2, , log t^i^^ + Y!<+uc^At log ? 

n i2 ^0 i ^0 

+E'el+i,a.4 log tl— ^ + e^t4 log ^ 

i CO ?0 

+E'^'('^*+1 ~ ^i) + ^"^('^'fc - + ^m'(T - Tk-l) + hai/3 - Tn) + /lajTl (Cll) 

i 

Here the double primed summation extends over the neighbors of the and r, which includes 

7'n,n,Tfc_i,Tfc. 

We are now ready to perform the RG calculation, by integrating out degrees of freedom 
in the range (C; C + ^0- This again leads to renormalization contributions of three kinds: a) 
simple rescaling; b) integrating out non-neutral pairs; and c) integrating neutral pairs. The 
contributions from simple rescaling of the cutoff, as well as from integrating out the primed 
kink-pairs, give rise to scaling of the fugacities, stiffness constants and the fields, which are 
the same as those in the partition function. 
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On the other hand, integrating out the pairs formed between a doubly primed kink 
and its neighbor in the primed kinks will contribute directly to the renormalization of the 
correlation functions. For the Green's functions, there are four terms contributing to the 
renormalization of the propagators. Contributions from the pair of kinks before and after 
the external time r cancel unless the original states before and after r are different. 

Collecting all these contributions, we arrive at the results 

d{addln{T / ^)) _ j^p^p^ „2 v„I 

d{acJn{T/^)) _ 2 _ m 2 \j 

u,i,ibi^ m A 7 ^ 

d{aJn{T/^)) _ 

where add, cxcc-, ctdc are the exponents for Gdd,Gcc,Gdc respectively, in the algebraic depen- 
dence on imaginary time. Integrating from ^ = .^o and C, = t gives rise to the correction to 
the Green's function due to the fugacities. 

The two particle correlation functions can be calculated following the same procedure. 
The results are given in section V A. 
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FIGURES 

FIG. 1. Hopping sequences in the atomic representation. Here r,, for i = l,...,n, labels the 
imaginary time at which the atomic configurations hops from \ai > to |aj+i >. 

FIG. 2. Phase diagram of the generahzed Anderson model in the mixed valence regime. Here 
€t and €j label the charge and spin stiffness constants defined in the text. The vertical thick line, 
the horizontal thick line, and the dashed line are the boundaries between the strong couping, weak 
coupling, and intermediate coupling mixed valence states. 

FIG. 3. Zero temperature transition sequences when the interactions are varied such that within 
the et — ej plane of Fig. |2|, (a) a one stage transition takes place from the weak coupling to the 
strong coupling state; and (b) two stage transitions take place from the weak coupling, through 
the intermediate coupling, into the strong coupling state. Here Ec, ed and ec2 label points on the 
vertical thick line, the horizontal thick line and the dashed line respectively in Fig. |2[ 

FIG. 4. Crossover diagram in terms of temperature (T) and chemical potential (/i) (a) for the 
strong coupling case. The dashed lines correspond to the crossovers, and A* is the renormalized 
resonance width; (b) for both the weak coupling and the intermediate coupling cases. Here Hc 
labels the critical point. 

FIG. 5. Electron density (n) as a function of chemical potential {n) (a) for the strong coupling 
case and (b) for both the weak coupling and the intermediate coupling cases. 
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FIG. 6. Schematic interaction— density phase diagram for two ways of varying the interactions. 
Here ei and €2 are the same as those in Fig. (a) For ei above Cc, non- Fermi hquid phase 
with both charge and spin excitations incoherent occurs tih zero temperature within the sohd-hne 
shaded density range, and above a crossover temperature within the dashed-hne shaded density 
range. Fermi hquid phase occurs otherwise, (b) For €2 above ed, non- Fermi hquid phase of (a) 
occurs. For €2 between ed and ec2, non- Fermi hquid phase with incoherent charge excitations but 
coherent spin excitations occurs tiU zero temperature within the sohd-vertical-hne shaded density 
range, and above a crossover temperature within dashed-vertical-hne shaded density range. Fermi 
hquid phase occurs otherwise. 
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